Majority-vote model on (3,4,6,4) and (3 4 ,6) Archimedean lattices 
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On Archimedean lattices, the Ising model exhibits spontaneous ordering. Two examples of these 
lattices of the majority-vote model with noise are considered and studied through extensive Monte 
Carlo simulations. The order/disorder phase transition is observed in this system. The calculated 
, values of the critical noise parameter are q c = 0.091(2) and q c = 0.134(3) for (3,4,6,4) and (3 4 ,6) 

Archimedean lattices, respectively. The critical exponents f3/v, f/v and 1/v for this model are 
0.103(6), 1.596(54), 0.872(85) for (3,4,6,4) and 0.114(3), 1.632(35), 0.98(10) for (3 4 ,6) Archimedean 
lattices. These results differs from the usual Ising model results and the majority-vote model on 
so-far studied regular lattices or complex networks. The effective dimensionalities of the system 
[D eS (3, 4, 6, 4) = 1.802(55) and D eff (3 4 ,6) = 1.860(34)] for these networks are reasonably close to 
the embedding dimension two. 
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I. INTRODUCTION 



■ 

The majority- vote model (MVM) |]J denned on regular lattices shows second-order phase transition with critical 
exponents (3, 7, v — which characterize the system in the vicinity of the phase transition — identical OlSISl with 
O \ those of equilibrium Ising model 0, |f| . 

On the other hand MVM on the complex networks exhibit different behavior 0, Gi H, El • Campos et al. investigated 
(*C) \ MVM on small- world network This network was constructed using the square lattice (SL) by the rewiring 
■ procedure. Campos et al. found that the critical exponents 7/1/ and (3/v arc different from these of the Ising model 
|{| and depend on the rewiring probability. Pereira et al. []| studied MVM on Erdos-Renyi's (ER) classical random 
graphs ^3 1 and Lima et al. Q also studied this model on random Voronoy-Delaunay lattice with periodic 
boundary conditions. Very recently Lima [ij studied the MVM on directed Albert-Barabasi (AB) network |TlT | and 
>^J , contrary to the Ising model m on these networks [T3I , the order / disorder phase transition is observed in this system. 



The calculated [3/v and j/v exponents are different from those for the Ising model || and depend on the mean value 
f^i ] of connectivity z of AB network. The latter was observed also for ER random graph Q. 
"r-7 ■ The results obtained by these authors 0, 0, H, 13 show that the MVM on various complex topologies belongs to 
different universality classes. Moreover, contrary to MVM on regular lattices Q, 0, the obtained critical exponents 
are different from those of the universality class to which the equilibrium Ising model belongs . 

In this paper we study the MVM on two Archimedean lattices (AL), namely on (3, 4, 6, 4) and (3 4 , 6). The topologies 
3?! of (3,4,6,4) and (3 4 , 6) AL are presented in Fig. The AL are vertex transitive graphs that can be embedded in a 
g ■ plane such that every face is a regular polygon. Kepler showed that there are exactly eleven such graphs. The AL 
q I are labeled according to the sizes of faces incident to a given vertex. The face sizes are sorted, starting from the face 
• • ■ for which the list is the smallest in lexicographical order. In this way, the square lattice gets the name (4,4,4,4), 
abbreviated to (4 4 ), honeycomb is called (6 3 ) and Kagome is (3,6,3,6). Critical properties of these lattices were 
investigated in terms of site percolation 0| in Ref. Top olog ies of all eleven AL are given there as well. Very 

^ ' recently, the Ising model on those AL was investigated in Ref. |16l. 

Our main goal is to check the hypothesis of Grinstein et al. — i-e., that non-equilibrium stochastic spin systems 
with up-down symmetry fall in the universality class of the equilibrium Ising model — for systems in-between ordinary, 
regular lattices (like square one H) and complex spin systems (like decorated with spin variables classical random 



graphs [3 or scale- free networks |9|) 
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FIG. 1: (3,4,6,4) (left) and (3 4 ,6) (right) AL. 
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FIG. 2: The magnetization M as a function of the noise parameter q, for L = 8, 16, 32, 64 and 100 sites for (a) (3, 4, 6, 4) and 
(b) (3 4 ,6) AL. 

With extensive Monte Carlo simulation we show that MVM on (3, 4, 6, 4) and (3 4 , 6) regular AL exhibits second-order 
phase transition with effective dimensionality D c ff « 1.8 and has critical exponents that do not fall into universality 
class of the equilibrium Ising model (although for (3 4 ,6) lattice they are quite close, yet lower). 



II. MODEL AND SIMULATION 



We consider the MVM pj defined by a set of "voters" or spin variables a taking the values +1 or —1, situated on 
every node of the (3,4,6,4) and (3 4 ,6) AL with N = 6L 2 sites. The evolution is governed by single spin-flip like 
dynamics with a probability Wi of i-th spin flip given by 



1 - (1 - 2q)Oi ■ sign ( ^ <jj 



(1) 



where the sum runs over the number z (z(3, 4, 6, 4) — 4 or z(3 4 , 6) = 5) of nearest neighbors of i-th spin. The control 
parameter < q < 1 plays the role of the temperature in equilibrium systems and measures the probability of parallel 
aligning to the majority of neighbors. It means, that given spin i adopts the majority sign of its neighbors with 
probability q and the minority sign with probability (1 — q) 0,Q,H,I3- 



To study the critical behavior of the model we define the variable m = J^iLi ■ In particular, we are interested 
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in the magnetization M, susceptibility x & n d the reduced fourth-order cumulant U 

M(q) = (\m\), (2a) 

X (q) = N «m 2 ) - (m) 2 ) , (2b) 

17(g) = 1 - <m 4 }/ (3(m 2 ) 2 ) , (2c) 

where (■ • • ) stands for a thermodynamics average. The results arc averaged over the iV run independent simulations. 
These quantities are functions of the noise parameter q and obey the finite-size scaling relations 

M = L-P/»f m {x), (3a) 

X = L^f x (x), (3b) 



d ^=L^f u {x), (3c) 

where v, (3, and 7 are the usual critical exponents, fi(x) are the finite size scaling functions with 

x = (q - q c )L x '" (3d) 

being the scaling variable. Therefore, from the size dependence of M and x we obtained the exponents (3/v and 7/1/, 
respectively. The maximum value of susceptibility also scales as L 1 ^. Moreover, the value of q* for which x has a 
maximum is expected to scale with the system size as 

q* = q c + bL- llv with b « 1. (4) 

Therefore, the relations (|3*c|) and Q may be used to get the exponent 1/v. We evaluate also the effective dimensionality, 
D c ff, from the hyperscaling hypothesis 

2/9/iz + 7/f = D cff . (5) 

We performed Monte Carlo simulation on the (3, 4, 6, 4) and (3 4 , 6) AL with various systems of size N = 384, 1536, 
6144, 24576 and 60000. It takes 2 x 10 5 Monte Carlo steps (MCS) to make the system reach the steady state, and 
then the time averages are estimated over the next 2 x 10 5 MCS. One MCS is accomplished after all the N spins are 
investigated whether they flip or not. We carried out N run = 20 to 50 independent simulation runs for each lattice 
and for given set of parameters (q,N). 



III. RESULTS AND DISCUSSION 



In Fig. |2]we show the dependence of the magnetization M on the noise parameter q, obtained from simulations on 
(3,4,6,4) and (3 4 ,6) AL with N ranging from 384 to 60000 sites. The shape of M(q) curve, for a given value of N, 
suggests the presents of the second-order phase transition in the system. The phase transition occurs at the value of 
the critical noise parameter q c . 

In Fig. |2| the corresponding behavior of the susceptibility x is presented. In Fig. 0]we plot the Binder's fourth- 
order cumulant U for different values of the system size N. The critical noise parameter q c is estimated as the point 
where the curves for different system sizes TV intercept each other [l^. ^From Fig. 0]we obtain q c = 0.091(2) and 
q c = 0.134(3) for (3,4,6,4) and (3 4 ,6) AL, respectively. 

In Fig. [S] we plot the dependence of the magnetization M* = M(q c ) vs. the linear system size L. The slopes of 
curves correspond to the exponent ratio j3/v according to Eq. (|3a|) . The obtained exponents are (3/v = 0.103(6) 
and 0.114(3) respectively for (3,4,6,4) and (3 4 ,6) AL. The exponents ratio 7/t/ arc obtained from the slopes of the 
straight lines with 7/2/ = 1.596(54) for (3,4,6,4) and 7/2/ = 1.632(35) for (3 4 ,6), as presented in Fig. 



4 



(3,4,6,4) 



(3 4 ,6) 



5 - 



?3 



2 - 









+ L=8 (a) 


- 


A 

□ 6 




x L=16 
o L=32 
□ L=64 


A 




A 
□ 


a L=100 


□ 
O 

□ 


o 

o 


o * 

o 


□ 

° 6 


ft 

o 

x 

X 


X 




* * x 6 e 


a * 

A x 






+ + + + + + + 


O + 

a + 

x + + 


+ + 


+ + 





6.0 



5.0 



4.0 



c 3.0 



2.0 



1.0 



0.0 





0.085 0.088 0.091 0.094 0.097 0.100 0.103 0.1 

q 











A 


+ L=8 

x I -1R 
x l_— I O 




A A 
□ 


o L=32 




□ 


□ L=64 






a L=100 




□ o o A 






A O O a 

OA 






o 5 


□ 

A 







6 □ 




g 

O x x 


x « □ 

x 






x □ 6 




o * 


X □ . 

*°„ 




o * 


+ + + + + + + + + + + + + «. 




Ax , , + 
OX + + + + 




x + + 

.8 , + + 





25 0.130 0.135 0.140 0.145 

q 



0.150 0.155 



FIG. 3: Susceptibility versus q for (a) (3,4,6,4) and (b) (3 4 , 6) AL 
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FIG. 4: The reduced Binder's fourth-order cumulant U as a function of q for (a) (3,4,6,4) and (b) (3 4 ,6) AL. 



To obtain the critical exponent l/i/, we used the scaling relation Q. The calculated values of the exponents 1/v 
are 0.872(85) for (3,4,6,4) (circle) and 1/v = 0.98(10) for (3 4 ,6) (square) (see Fig. 0). Eq. © yields effective 
dimensionality of systems D oS = 1.802(55) for (3,4,6,4) and D eS = 1.860(34) for (3 4 ,6). The MVM on those two 
AL has the effective dimensionality close to two contrary to ER classical random graphs (0.99 < D e ff < 1.02) or 
directed AB networks (0.998 < D cff < 1.018) with roughly the same nodes connectivity [z = 4) as for (3,4,6,4) 
and (34,6) AL. 

The results of simulations together with data for ER and AB networks with mean connectivity z = 4 are collected 
in Tab. [Q 
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FIG. 6: Critical behavior of the susceptibility x(N) at q — q c (circle) and q = q* (diamond) for (a) (3,4,6,4) and (b) (3 4 ,6) 
AL. 



IV. CONCLUSION 



We presented a very simple non-e quilibrium MVM on (3, 4, 6, 4) and (3 4 , 6) AL. On these lattices, the MVM shows 
a second-order phase transition. Our Monte Carlo simulations demonstrate that the effective dimensionality D e g is 
close to two, i.e. that hyperscaling may be valid. 

Finally, we remark that the critical exponents j/u, (3/v and \jv for MVM on regular (3,4,6,4) AL are different 
from the Ising model p| and differ from those for so-far studied regular lattices QJ3and for the ER classical random 
graphs and for the directed AB network However, in the latter cases |9( the scaling relations (j3J) must 
involve the number of sites N instead of linear system size L as these networks in natural way do not posses such 
characteristic which allow for N oc L d (d G Z) dependence |20|. For (3 4 , 6) AL, the critical exponents are much closer 
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FIG. 7: The exponents 1/v obtained from the relation for (3,4,6,4) (circle) and (3 4 ,6) (square) AL. 

TABLE I: Critical parameter, exponents and effective dimension for MVM model on (3, 4, 6, 4) and (3 4 , 6) AL. For completeness 
we cite data for scale-free AB and random ER graphs with an average connectivity 2 = 4. 



sl m 



AB [9] 



ER [71 



(3,4,6,4) 



(3 4 ,6) 



Qc 

l/u 
D eB c 



0.075(10) 
0.125(5) 
1.73 
1.70 

1/0.99(5) 



0.431(3) 
0.447(1) 
0.104(2) 
0.888(9) 

0.998(3) 



0.181(1) 

0.242(6) 

0.54(1) 

0.515(6) 

0.59(7) 

1.02(2) 



0.091(2) 

0.103(6) 

1.596(54) 

1.663(32) 

0.872(85) 

1.802(55) 



0.134(3) 

0.114(3) 

1.632(35) 

1.681(10) 

0.98(10) 

1.860(34) 



"obtained using x(N) at q = q c 
''obtained using x{N) at q = q* 

"obtained using ratio ■y/u given by dependence x(^0 at q = q c 



to those known analytically for square lattice Ising model, i.e. (3=1/8 = 0.125, 7 = 7/4 = 1.75 and v = 1, but except 
for v they differ for more than three numerically estimated uncertainties. As those uncertinities were estimated basing 
only on statistical analysis we cannot excluding that critical exponents may reach SL Ising exponents asymptotically 
in the thermodynamic limit. 
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